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Current theoretical results for the three-dimensional Navier-Stokes equations only 
guarantee that solutions remain regular for all time when the initial enstrophy 
(||Z)uo|| 2 : = / |curlwo| 2 ) is sufficiently small, ||-Dwo|| 2 < Xo- In fact, this smallncss 
condition is such that the enstrophy is always non-increasing. In this paper we pro- 
vide a numerical procedure that will verify regularity of solutions for any bounded 
set of initial conditions, ||-Dw,o|| 2 < Xi- Under the assumption that the equations are 
in fact regular we show that this procedure can be guaranteed to terminate after a 
finite time. 
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1. Introduction 

The global regularity of solutions of the three-dimensional Navier-Stokes equations 
is one of the most intriguing open problems in mathematics. Despite the efforts 
of many eminent mathematicians over the last 150 years, the question whether 
the most fundamental model in hydrodynamics is or is not self consistent remains 
unanswered. Even when the flow is not disturbed by any external forces and the 
problem of boundaries is avoided by considering periodic boundary conditions (or 
the whole of K 3 ) the problem is unresolved. 

More precisely, consider the Navier-Stokes equations for the flow of an incom- 
pressible fluid 

du 

— -uAu+Ou- V)u + Vp = V-u = u(0)=u . (1.1) 
at 

with periodic boundary conditions on Q — [0,-L] 3 , where v > is the kinematic 
viscosity. It is known that weak solutions of these equations exist that are defined 
for all t > 0. However, we do not know whether these solutions are unique and 
regular. On the other hand, we know that if the initial enstrophy is finite then the 
solutions are regular for t sufficiently small (t < T±) and sufficiently large (f > T2), 
where T\ and T2 depend respectively on the enstrophy and kinetic energy of initial 
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conditions: 

Ti ~ ||L»uoir 4 and T 2 - ||u || 2 , (1.2) 

where 

\\ u o\\ 2 = / |w(a;)| 2 dx and ||Du || 2 = / |curlu(x)| 2 da;. 
Jo Jn 

It is an open problem whether or not weak solutions develop any singularity at 
some time T\ < t < T 2 . 

In the light of (1.2) it is natural to conjecture that some restrictions on en- 
strophy or energy of initial conditions might provide us with regularity of weak 
solutions. However, nothing is known about regularity of weak solutions only under 
the assumption that energy of the flow is small. In other words, there is no theoret- 
ical evidence that there exists an e > such that all initial conditions with kinetic 
energy 

l|w<o|| 2 < £ 

do not blow up. On the other hand, under the stronger assumption that the enstro- 
phy of the initial condition uo is sufficiently small, 

\\Du \\ 2 <ei 

one can prove that uo gives rise to a strong global solution. The value of E\ for which 
the argument holds is very small (see Section 2). This is due to the fact that the 
proof of regularity of solutions arising from small initial data (sketched in the next 
section) actually rules out all initial conditions for which the enstrophy is increasing 
at the beginning of the flow. Thus it is desirable to extend this theoretical result to 
a larger class of initial conditions that allow for more complicated and physically 
realistic behaviour of the flow. 

The main aim of this paper is to present a numerical method which may provide 
us with such an extension. We show that, under the assumption that the equations 
are well-posed, it is possible to verify numerically the regularity of solutions for all 
initial conditions with ||-Dwo|| < R for any choice of R > 0. In Section 2 we set 
up the problem more formally, and recall some classical results on the regularity of 
solutions. Section 3 reproduces results from Dashti & Robinson (2006) that can be 
used to verify that a single initial condition u 6 V gives rise to a strong solution 
of some bounded interval [0, T]. Section 4.2 shows that for initial conditions in 
any bounded ball in V it is sufficient to guarantee that solutions are regular on 
some finite time interval [0, T], since all solutions are regular after some finite time. 
Section 5 is the core of the paper, and shows that one can verify that all initial 
conditions in some bounded ball in H 2 give rise to regular solutions - this requires 
no additional assumptions. The final section shows his this can be used to verify 
the same property for all initial conditions in a bounded ball in H 1 , although we 
note here that this requires the assumption of regularity for all initial conditions in 
H 1 . 



2. Setting of the problem 

As enstrophy and kinetic energy play key role in the problem of existence and 
regularity of solutions it is worth recalling two basic facts concerning them. First 
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of all, it follows from the energy inequality 

^H 2 + I /p?«|| 2 <o (2.i) 

that the kinetic energy of the flow decreases to zero as time goes to infinity (the 
flow is dissipative). Moreover, we have 



Jo 



T 

2j„ / ..-l|L./n\l|2 



\\Du(s)rds<^\\u(0)r (2.2) 



which shows that mean value of ||-Du|| 2 over the interval [0, T] is bounded in terms 
of the initial kinetic energy and decreasing at least as fast as T _1 . 

Throughout the paper we use standard notation for the Lebesgue and Sobolev 
spaces L 2 , H 1 and H 2 . The set of all smooth, periodic and divergent free functions 
on Q, with u{x) Ax = we denote by V. The closure of V in the L 2 and H 1 norms 
we denote by H and V, respectively. The eigenfunctions of the Stokes operator A 
defined by 

A{u) = -PA(u), 

where P is the Leray projector (the orthogonal projection of the space L 2 onto H), 
we denote by Wk and their corresponding eigenvalues by It follows from the 
theory of compact self-adjoint operators in Hilbert spaces that we can order the 
eigenvalues so that 

< Ai < A 2 < A 3 < ... — > oo. 

Since we consider periodic boundary conditions we can in fact specify the values of 
Afe exactly, and in particular we have 

4^ 2 

~~ T 2 ' 

In what follows we will concentrate on a case L = 2ir for which Ai = 1. This is 
convenient in the light of the following Poincare inequalities 



it 



2 



<X^\\Du\\ 2 , \\Du\\ 2 < X^\\Au\\' 



which we use frequently in what follows. 

As mentioned above, solutions to (1.1) with sufficiently small initial data are 
regular for all time. With zero forcing the condition guaranteeing regularity is given 
by following theorem. 

Theorem 2.1. If the initial condition u n for problem (1.1) satisfies 

\\Du Q \\ 2 < c-VVAi 72 , (2.3) 
where c is some absolute constant^, then ua gives rise to a global strong solution. 



f For later use it is worthwhile to note that c = 27fc 4 /16, where k is the constant in the 
inequality 

< kWDuUDv^^WAvf^lAwl 



\ I [(« ■ V)u] • Awdx 
\Jn 



In the cubic geometry here we can take k = 9 x 2 15 / 4 , see Dashti & Robinson (2006). 



Article submitted to Royal Society 



4 



J.C. Robinson & W. Sadowski 



The proof is based on the inequality 

^WDufK^WDut-yMDuf (2.4) 

which may be obtained from the enstrophy inequality with some help from the 
Sobolev embedding theorem. What is crucial is that for initial conditions satisfying 
(2.3) the right hand side of (2.4) is non-positive. In consequence ||-Du|| 2 will never 
exceed the value c~ 1//2 ^ 2 aJ / ' 2 . (For details of the proof see, for example, Theorem 
* in Constantin & Foias (1988).) 

Below we treat the question of regularity for periodic boundary conditions on 
[0, 2ir] 3 and for viscosity v = 1. We now show, using a simple scaling argument, that 
a proof of regularity for this particular case would imply regularity for all [0, L} 3 
and all v > 0. Indeed, starting with the equation 

— -vAu+(u-V)u + Vp = V-m = 0, ie[0,L] 3 , 

we set 

u = Lu/2-kv, x — 2ttx/L, t = 4tt 2 vt / L 2 , and p — L 2 p / 'An 2 v 2 , 
and then obtain 

du 

- A z u + {u ■ V z )u + Vsp = (2.5) 

at 

for x e [0,2tt] 3 . 

For this non-dimcnsionalizcd form of the Navier-Stokes equations we can get 
an explicit value Ry for the radius of the ball in V in which standard results 
guarantee that every initial condition gives rise to a global strong solution. Indeed, 
since c = 27(9 x 2 15 / 4 ) 4 /16 we get 

Rv - c- 1 / 4 ^ 1 / 2 = c-i = 2/(9 x 2 15 / 4 x 27 1/4 ) w 0.007. 

Nothing is known about regularity of solutions starting from initial conditions lying 
in a ball in V with R > Ry, and so the theoretical results currently available 
really are about "small data" . In this paper we present a method of numerical 
verification of regularity for more general, physically reasonable initial conditions. 
More precisely, denote by SSy[R) the ball of radius R in V. 

We prove below under the assumption of regularity that for any fixed value of 
R > the following statement can be verified numerically in finite time. 

Statement 2.2. Every initial condition uq € 3§y{R) gives rise to a global strong 
solution of problem (1.1). 

Of course, regularity of solutions with initial conditions in some fixed ball in 
V for the non-dimensionalized form of the Navier-Stokes equations implies parallel 
results for other cubic domains and different values of the viscosity. Indeed, scaling 
initial condition changes its enstrophy as we have 

du L 2 du 
dx A-k 2 v dx 
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So 




and since Ai = we see that 

im 2 = ^im 2 

which is in agreement with theorem 2.1. 

Thus a complete result on the regularity of solutions of the non-dimensionalizcd 
Navier-Stokes equations would immediately imply the same result for any L and 
any v, while regularity only for ||-Du|| 2 < R in the non-dimensionalized equation 
would provide regularity for all ||-Du|| 2 < 2ttRv 2 /L. 



3. Numerical verification of regularity for a single u e V 

Before presenting a method of verifying regularity for all initial conditions in bounded 
subsets of H 2 (and then H 1 ) we need a result that allows us to check numerically 
the regularity of a particular solution arising from some given initial condition. In 
other words, we have to be able to make a judgement based only on numerical cal- 
culations whether a given initial condition uq gives rise to a regular solution u. The 
first step to this end is a standard one: We construct the Galerkin approximations 
u n of u by solving the equation 

— + Au n + P n B(u n ,u n ) = 

with u n (0) = P n uo, where P n is orthogonal projection on the space spanned by 
first n eigenfunctions of the Stokes operator. 

We can deduce solely from the properties of these Galerkin approximations 
whether they approximate a regular solution on a finite time interval. A method 
for this was provided by Chernyshenko et al. (2006) under the assumption that 
ito 6 H m with m > 3. Dashti & Robinson (2006) adapt these arguments to show 
that the same is possible given a regular solution arising from an initial condition 
in V, and wc now briefly outline the main idea of this method. 

The first ingredient is a robustness result for strong solutions: 

Theorem 3.1. There exists a constant c > such that if u e V gives rise to a 
strong solution u(t) of (2.5) on [0,T*] then so does any initial condition v and 
forcing g{t) with 

\\D(v -u )\\+£ \\Dg(s)\\ ds < C (T*)- 1/4 cx P {-c £ \\Du{s)t + ||£>«( a )|| \\Au{s)\\ As 

(3.1) 

(The constant c is related to the constants arising in certain Sobolev embedding 
results.) 

This has the following corollary, key to all that follows: 
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Corollary 3.2. Suppose that v e L°°(0,T*;T/) fl L 2 (0, T*;D(A)) is a numerical 
approximation of u such that 

(ill 

— + Av + B(v, v) e L^O, T*; V) n L 2 (0, T*\H) 



and 

\Dv(0)-Du \ 



o 



— (s) + A«(s) +£(«(«),«(«)) 



< c(T*)- 1/4 cxp i^-c£ \Dv(s)\ 4 + \Dv{s)\\Av(s)\ds^J . (3.2) 

Then u is a regular solution of (2.5) with u £ L°°(0,T; V) fl £ 2 (0,T; D(A)). 

Thus to verify that Mo 6 V gives rise to a strong solution on [0, T*], it suffices 
to show that for sufficiently large n the Galerkin solution u n will satisfy (3.2) in 
Corollary 3.2. If this is the case then on can compute ever more accurate Galerkin 
approximations until (3.2) is eventually satisfied. If the approximated solution is 
regular then this method has to be successful since it is the case that the Galerkin 
approximations of regular solutions must converge in a sufficiently strong sense: 

Theorem 3.3. If uq E V gives rise to a strong solution u{t) on [0, T*] then the 
Galerkin approximations u n converge to u strongly in L°°(0, T; V)CiL 2 (0, T; D(A)). 

Of course, this neglects two issues: first, that the numerical solution will not be 
a continuous function of t, but given at discrete time points (cf. comments in 
Chernyshenko et al. (2006)). However, since Corollary 3.2 says nothing of the prove- 
nance of v, one can construct v via linear interpolation between these time points. 
A further complication is that the integrals in (3.2) cannot be computed exactly 
and have to be approximated; but this is not a major obstruction. 



ds 



4. Finite time of calculation 

Our goal in the following two sections is to prove, under the assumption of regularity, 
that the following statement may be verified numerically in a finite time. 

Statement 4.1. Consider the non-dimensional 3d Navier-Stokes equations on 
[0, 27r] 3 with zero forcing, 

— - vAu+ (u ■ V)w + Vp = V-w = u(0)=u . (4.1) 

Then every initial condition in a ball in V 

B V {R) = {u : ||D«o|| < R} 

gives rise to a strong solution that exists for all t > 0. 

Since the method in the previous section only works to verify the regularity of 
a given solution on some finite time interval, it is vital to show that this is all that 
is required. 
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Lemma 4.2. Any (weak) solution of (4-1) becomes regular after a finite time 

T* = c 1 / 2 \\u \\ 2 <c 1 / 2 \\Du \\ 2 , 
where c is the absolute constant from theorem 2.1. 
Proof. Since 

/ pn( S )|| 2 d S <||n || 2 , 
Jo 

it follows that on any time interval [0, T] with T > T* as above, there must exist a 
t such that 

\\Du(t )\\ 2 <c-^ 2 , 

and hence from theorem 2.1 it follows that the solution is regular for alH > T* . □ 

Thus if we are interested in regularity of a particular ball 3§v{R) we have to 
check regularity of solutions only on some fixed interval of time: now we can proceed 
with the proof of the main theorem of this paper. 

The most straightforward approach would be to verify the regularity of solutions 
arising from some finite sequence of initial conditions u^ n , where 1 < n < N 
and iV is large enough. Then Corollary 3.2 would imply regularity of all initial 
conditions in balls with radii Si, 82, Sn centred at uqi, uq2, uqn respectively. 
If these balls cover all of 33r then the Statement 4.1 would be verified. However, 
several obstacles stand in the way of such a happy ending. First of all, a ball in 
infinite dimensional space is never compact and one cannot cover it with any finite 
number of smaller balls. Even if we were interested only in checking regularity in 
some finite-dimensional subspace of V we would still encounter a serious problem: 
without a proof of some general uniform lower bound on the sequence Sk we might 
not be able to construct a proper covering (e.g. it is not possible to cover the unit 
interval with a sequence of intervals of length S n = 3~ n ). 

To circumvent these difficulties we will prove in the next section an auxiliary 
result under a stronger assumption on initial condition. 

5. Numerical verification of regularity for all initial 
conditions in a ball in H 2 

In this section we show that it is possible to verify numerically the following state- 
ment, which requires more regularity of the initial condition. 

Statement 5.1. Consider the 3d Navier-Stokes equations on [0,27r] 3 with zero 
forcing, 

— - Au+ (u ■ V)u + Vp = V-u = u(0) = u o . (5.1) 
Then every initial condition in a ball 

&h*(S) = {uo:\Auo\<S} 
gives rise to a strong solution that exists for all t > 0. 
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Theorem 5.2. Statement 5.1 can be verified numerically in a finite time. 

This theorem is proved via a series of subsidiary results. From lemma 4.2 and 
the Poincare inequality it follows that in order to verify Statement 5.1 we need only 
prove regularity on the finite time interval [0,T*]. It is therefore very useful to have 
the following result, which gives uniform bounds over the H 1 and H 2 norms of any 
solution with m G £§h 2 {S)- 

Proposition 5.3. Suppose that Statement 5.1 is true. Then for any T > there 
exist constants Ds(T) and Es(T) such that for any solution u arising from an 
initial condition uo G £$h 2 (S) we have 

r T 

sup \\Du{s)\\<D s and / || Au(s)\\ 2 ds < E s . (5.2) 

0<t<T Jo 

Proof. This result follows from the proof, but not the statement of Theorem * in 
Foias & Temam (19**), see also Theorem 12.10 in Robinson (2001). Here we give 
a sketch of the relevant part of the argument. If the first bound in (5.2) does not 
hold then there must be sequences uo n with ||Auo n || < S and t n G [0, T] such that 
\\Du n (t n ) || — > oo as n — > oo. By extracting appropriate subsequences we can assume 
that t n — > T G [0, T] and that UQ. n — uq in D(A): in particular, ||Auo|| < S and 
uon —* Mo in H. As in the standard proof of the existence of weak solutions, one 
can show that the solutions u n (t) of 

du n /dt + Au n + B(u n , u n ) = u(0) = u Qn 

are uniformly bounded in L°°(0,T;H) n L 2 (0,T;V); taking the limit as n -> 
produces a weak solution v of 

dv/dt + Av + B(v,v) = with v(0)=v o . 

However, by assumption, this equation possesses a strong solution which is unique 
in the class of weak solutions, and so in fact v G L°°(0,T;V) n £ 2 (0, T; D{A)). 
One can now use this regularity of v along with standard estimates on the equation 
for the difference v — u n to show that u n — > v strongly in L 2 (0, T;V). Since L? 
convergence on an interval implies that there exists a subsequence converging almost 
everywhere, at each such point of convergence, s, we must have ||Dm„(s)|| < So := 
1 + |H|l°°(o.T;V) for n sufficiently large. Standard estimates can be used to show 
that there is a time r such that 

\\Du n (s + t)\\ < 2(1 + S ) for all 0<t<r 

for all such s and n. By covering [0, T] with a finite number of intervals [sj, Sj + r], 
where u n (sj) — ► u(sj) in V", one can deduce that in fact ||£)u„(f)|| < 2(1 + So) for 
all t G [0, T] and obtain a contradiction. The second bound in (5.2) then follows 
from the inequality 

^||D«|| 2 +||A«|| 2 < C ||£>«|| 6 
which can be obtained by standard methods. □ 
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Now we are ready to prove that if an initial condition uq G 3§h 2 (S) gives rise 
to a regular solution then there exist a ball in V of some radius 5 depending on 
S but not on uq such that all initial conditions in that ball give rise to regular 
solutions. In other words, we may say that if Statement 5.1 is true, then all points 
of &h 2 (S) are separated in V from any initial conditions giving rise to a solution 
with singularities (if such solutions exist) by at least some fixed distance 6. Most 
importantly for us it means that there is a uniform lower bound on the radius of 
the "balls of regularity" whose existence is guaranteed by theorem 3.1, centered at 
points of SS H i (S) . 

Proposition 5.4. Suppose that Statement 5.1 is true. Then there exists a 5 = S(S) 
such that if uq G SS h q. (S) then the value of 5 required by the regularity check of 
Corollary 3.2 is bounded below by some 8: any initial condition vq G V with 

||D«o-D«o|| < S 
also gives rise to a strong solution on [0, T*] . 

Proof. Given the result of Proposition 5.3, the integrand on the right-hand side of 
(3.1) is bounded above: 

J:= f T \\Du(s)\\*+\\Du( S )\\\\Au(s)\\ds<I s :=T*Dl ; + (T*D 2 s ) 1 /ZEs, 
Jo 

and so the right-hand side of (3.1) is bounded below, 

cexp(— cl) > 5 := ccxp(— els). 

It follows that there exist a lower bound on the radius of "ball of regularity" pro- 
vided by theorem 3.1 □ 

Now we use a compactness argument. Loosely speaking, we could cover SBjji {S) 
by open balls of radius 5 centred at every point of 3§h2 (S) and then choose finite 
sub-covering (in topology of V). This argument, based on the Rcllich-Kondrachov 
theorem on the compact embedding of H 2 in H 1 , is formulated more precisely in 
the following lemma which also gives the explicit form of points in SS^i (S) that we 
need to test for regularity numerically. 

Lemma 5.5. Given S > there exist N, M such that every v G 3§ H 2{S) can be 
approximated to within 5 (in the V norm) by elements of the set 

N 

Um.n = {uq : Mo = ^^ a j w j, with aj — aj/2 M , aj G Z; ||Au || < S} 
where {wj} are the Stokes eigenf unctions. 

Note that the proof shows that one can take N such that Ajv + i > 2S 2 /S and M 
such that 2- M < 5/2. 
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Proof. Take any v € 33 H i(S). Since \\Av\\ < S, we have 

oo 

\\D(P n v - v)\\ 2 = Yl X k\(v,w k )\ 2 

k=n+l 

oo 

< Kli E A2 fe i(^)i 2 

k=n+l 

< KliWM 2 

Thus we can choose an N that for any n > TV we have 

||I>(PnV- V )|| 2 <A^S 2 <! 

and we can cover the finite dimensional ball Pjv(^_y2 (5*)) with a finite number of 
balls of radius r < | with centres in Um,n, where 2~ M < (5/2. It follows that for 
any w £ ^_y2 (5*) we can find uq £ Wm,jv such that 

||D«b - Du a \\ < \\D(v - P N v )\\ + \\P N v a - Du \\ < S. 

□ 

Finally we assemble all these ingredients to provide a proof of theorem 5.2. 

Proof of thereom 5.2. Given S > 0, let 5 be the lower bound on the radius of 
the "ball of regularity" provided by proposition 5.4. Then to prove regularity for 
all initial conditions in 3§h 2 (S) it suffices to check that the finite set of initial 
conditions in Um,n (the set defined in lemma 5.5) all give rise to strong solutions 
on [0, T*]. This can be done for each initial condition in turn using the procedure 
outlined in Section 3. □ 

Observe that the proof is based only on an assumption of regularity of solutions 
arising from initial conditions in .^§h 2 (S). Whether or not singularities may occur 
in solutions arising from initial conditions mo with ||Auo|| > S has no effect on the 
possibility of numerical verification of regularity in the ball £$ H 2 (S) . 

6. Numerical verification of regularity for all w £ ^y(-R) 

Now we will use the results of previous section to prove the main theorem of this 
paper 

Theorem 6.1. If the Navier-Stokes equations are regular, then statement 4-1 can 
be verified numerically in a finite time. 

Proof. We find a Tr > such that every solution with initial condition H-Dwoll < R 
is regular on the interval [0,Tr] and has an explicit bound on ||Am(Tr)||, using a 
result due to Foias & Temam (1988) on the Gevrey regularity of solutions. 
More explicitly, if we define 

00 , n 

e tA^ 2 = vlw 2 



= E h An/2 ^ 



n 

n=0 
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and set 

T=-L(l+||^ || 2 )- 2 , 

where K\ < 3266 is an absolute constant depending on certain algebraic and 
Sobolev embedding inequalities, Foias & Temam (1988) showed that for < t < 2r, 
one has 

\\A^e tAl/2 u(t)\\ 2 <2(l + \\Du r). 

By expanding A 1 / 2 e tA ' u in terms of the eigenfunctions of A (cf. proof of 
Lemma 2 in Friz & Robinson, 2001) one can show that 

^f\\Auf < \\A^ Al/2 u\\, 

and hence that 

||A U (r)|| 2 <X 2 (l + !|^ || 2 ) 5 . 

So it suffices to check for regularity of initial conditions in some ball of radius 
S := Ki(l + i? 2 ) 5/2 in H 2 . Indeed, if any initial condition u E & V {R) led to 
singularity so would some initial condition uq = u(t) e £$h 2 (S). So the regularity of 
solutions arising from £$ H i{S) - which according to 5.2 can be verified numerically 
- implies the regularity of solutions arising from {%y (R) . □ 

Observe that in order to prove Statement 4.1 we need to assume that all solutions 
of the Navier-Stokes equations are regular. Indeed, we cannot exclude the possibility 
that some singularity may occur in a solution arising from an initial condition uq 
with ||-Dwo|| > R but for which ||^4u(r)|| < S. Then, of course, we would not be 
able to "verify" numerically a false statement about regularity of solutions arising 
from a ball 33 H i (S) in spite of the fact that solutions arising from a ball SSy (R) 
would be regular. 

7. Conclusion 

We have shown that if the Navier-Stokes equations are regular on a cubic periodic 
domain, then regularity for any bounded set of initial conditions in V (bounded 
enstrophy) can be verified numerically using an algorithm that is guaranteed to 
terminate. While this falls far short of a numerical proof of regularity, it is at least 
theoretically possible to check the validity of the equations for all initial conditions 
that one could hope to realise in experiment. 

Of course, in this approach to regularity the outstanding open problem is to 
show that in fact regularity for all initial conditions in some bounded ball 5§y(R) 
is in fact sufficient for regularity for all initial conditions, but this seems to be 
an extremely strong statement and it appears unlikely that this could be proved 
analytically without leading to a full proof of regularity. 
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